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Abstract 

The effect of viscous drag and surface tension on the nonlinear two fluid interfacial struc- 
tures induced by Rayleigh -Taylor instability and Richtmyer-Meshkov instability are investi- 
gated. Viscosity and surface tension play important roles on the fluid instabilities. It is seen that 
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the magnitude of the suppression of the terminal growth rate of the tip of the bubble height 
depends only on the viscous coefficient of the upper (denser) fluid through which the bubble 
rises and surface tension of the interface. However, in regard to spike it is shown that in an 
inviscid fluid spike does not remain terminal but approaches like a free fall under gravity as the 
Atwood number (A) increases. In this respect there exits qualitative agreement of our results 
with simulation result as also with some earlier theoretical results. Viscosity reduces the free fall 
velocity appreciably and it becomes terminal with increasing viscosity. Results obtained from 
numerical integration of the relevant nonlinear equations describing the temporal development 
of the spike support the foregoing observations. 
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I INTRODUCTION 

There are many causes for instability of the interface between two fluids. Under the gravitational 
force when a denser fluid overlies a lighter fluid, the instability occurs and it is called Rayleigh- Taylor 
instability(RTI). Ritchmeyer-Meshkov instability (RMI) is another type of instability which occurs 
whenever a shock front crosses the interface of two materials of different shock impendence (the shock 
must enter from the low independence interface side). Both these instabilities play important role in 
the ablation region at compression front during the process of inertial confinement fusion, supernova 
remnant formation or shock tube experiments in the laboratory [1]. In the nonlinear regime the fluid 
interface forms a finger shape structure. The structure is called a bubble (spike) if the lighter (denser) 
fluid penetrates into the denser (lighter) fluid. Under astrophysical conditions such structures may 
cover enormous range of spatial distribution. Examples are suggested to be provided by pillars 
("elephant trunk") of Eagle Nebula which is identified with the spike of a heavy fluid penetrating 
a lighter fluid[2-5]. Also sudden increase in the height of the ionospheric F 2 layer is caused by RTI 
mechanism as suggested by some observational data[6]. Layzer[7] was first to describe the bubble 
formation in a potential flow model. This model which is based on an approximate description of 
the flow near the bubble tip describes its nonlinear growth[8-13, 26-28]. 

The effect of viscosity on Rayleigh- Taylor instability and Richtmyer-Meshkov instability shows 
significant importance for increasing wave number k as uk 2 where v is the kinematic coefficient of 
viscosity. This effect is further enhanced as v increase with the temperature for a plasma [14]. The 
importance of this feature has been discussed for supernova remnant in Ref. [15]. In the domain of 
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linear theory, the effect of viscosity and surface tension on RTI was described in depth by Chan- 
drasekhar [16]. The same effect was described by Mikaelian[17-18] for RMI and for RTI in finite 
thickness [19]. Higher magnitude of the growth rate suppression due to viscosity was obtained in the 
weakly nonlinear theoretical study of Carles and Popinet[20]. Later, a time dependent expression of 
the reduced linear theory growth rate of RTI and RMI due to the combined effect of viscosity and 
binary mass diffusivity between the fluids was arrived at in linear theory by Robey[21]. The effect of 
surface tension under the weakly nonlinear theory was analyzed by Gamier et al. [22] and later Roy 
et al.[29] studied the same effect in nonlinear theory. 

The present paper reports the combined effect of viscosity and surface tension on the two fluid 
nonlinear interfacial finger like structures resulting due to RTI and RMI. We have analyzed the 
problem based on Layzer's approach. It is seen that, in absence of surface tension the lowering of 
the asymptotic velocity of the tip of the bubble which is formed when the lighter (the lower) fluid 
penetrates into the denser (the upper) fluid and thus encounters the viscous drag due to the denser 
fluid, depends only on the viscosity coefficient of denser fluid. However, in presence of surface tension, 
the asymptotic velocity of the tip of the bubble and nonlinear perturbed surface are oscillating under 
certain conditions. It has been shown that, for RTI this oscillation depends only on the surface 
tension but for RMI it depends on surface tension as well as viscosity. 

On the other hand it is shown that in an inviscid fluid the spikes do not remain terminal as 
obtained in theoretical works[10,13]. It is only discussed in ref.9 that the RTI spike was shown to 
have a free fall ~ |gt 2 where g is acceleration due to gravity and the RMI spike to have a constant 
velocity of fall but that too only for Atwood number A = 1. In the present analysis the free fall 
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behavior is seen to hold for all values of Atwood number A. The depth of the tip of the spike below 



A — > 1. Similar result is found to hold for RMI. The effect of viscosity is seen to reduce the spike 
velocity appreciably and as the viscosity coefficient increases the spike velocity tends to become 
approximately terminal. Such behavior of both RTI and RMI spike for for inviscid and viscous fluids 
was not found earlier. 

The paper is organized as follows: In Section II, we describe the basic fluid equation with the 
assumption that the motion is irrotational and fluid is incompressible. Also the kinematic and 
dynamical boundary condition are derived in Sec. II. The analytical expressions of the viscosity and 
surface tension induced asymptotic velocity of the tip of the bubble and associated numerical results 
are given in Section III. Section IV is devoted to temporal development of the spike. The results are 
summarized in Section V. 



II BASIC EQUATION AND BOUNDARY CONDITIONS 



The x — y plane (z = 0) is assumed to be the unperturbed interface between the denser fluid of 
density p h (region z > 0) and lighter fluid of density pi (region z < 0). The variables with subscript 
h and I represents denser and lighter fluid, respectively. Gravity g is taken to point along negative 
z-axis. After perturbation the finger shape interface is assumed to take up a parabolic form, given 



the unperturbed interfacial surface is ~ ^ 



2 where 7 is a dimensionless constant which — > 1 as 



by 



z = r)(x, t) = r) (t) + r] 2 (t)x 2 



(1) 
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For a bubble: i] > and r] 2 < 0, (2) 

For a spike: i]q < and r^ 2 > 0. (3) 

Following Goncharov [13], the velocity potentials describing the irrotational motion for the denser 
and lighter fluids are assumed to be given by 

<p h (x, z, t) = a x (t) cos {kx)e~ k{z ^ o{t)) ; z > (4) 

(j>i(x,z,t) = b {t)z + b^t) cos {kx)e kiz - r,o{t)) ; z<0 (5) 

where k is the perturbed wave number. 

To find the five unknown functions i.e. r) (t) ,r)i(t) ,ai(t) ,b (t) and bi(t), we require as many equa- 
tions obtained from the kinematical and dynamical boundary conditions describing the dynamics. 

We first turn to the kinematical boundary conditions corresponding to the interfacial surface 
perturbations represented by eq.(l): 

Vx(Vh)x ~ Vx(vi) x = (Vh)z ~ {vi) z (6) 
Vt + Vx(Vh)x = (Vh)z (7) 

Substituting eq.(6) and eq.(7)from eq.(l) and for (v h (i)) x = and ( v h(i))x = from 

eq.(4)and eq.(5) and expanding in powers of the transverse coordinate x neglecting terms 0(x l )(i > 
3), we obtain the following relations which are equivalent to the kinematic boundary conditions eq.(6) 
and eq.(7): 
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f^-i(6 6 + l)6 (9) 

bo = -^rka, (10) 

6 & - 2 

61 = ^7 r«i ( n ) 



where £1 = /cr^; 6 = ^2/^; 6 = & «i; T = t^kg (12) 

£1 and £ 2 are, respectively, the non-dimensionalized displacement and curvature at the tip of the 
nonlinear structure, ^/k is its velocity and r is the non-dimensionalized time. Eq.(8) and eq.(9) are 
the first two of the three time development equations needed to describe the time evaluation of the 
nonlinear structure (the other two viz bo(t) and b\{t) are provided by eq.(10) and eq.(ll)) . 

For the constant density fluid, the equation of continuity is V.v = 0, which becomes V 2 = for 
potential flow. So, for a fluid with uniform viscosity having coefficient of viscosity Hh{i), the viscous 
term drops out (fj,V 2 v = /zV(V 2 0) = 0) in 
dv 



P 



m + < lXV >" 



-Vp + /xW- pgz (13) 



and we arrive at the first integral of the momentum equation. 

- Ph(i)—^ + WMO) 2 + PHi)9 z = ~Ph(i) + fh(i)(t) (14) 

The net stress [16] at two fluid interface including that due to viscosity is 

^(«) = -p fc (0 + 2^(0^^ (15) 
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Plugging the dynamical boundary condition P h — P h = —T/R at the interface z(x, t) = rj(x : t), where 
R is the radius of curvature and T is the surface tension of the perturbed interface, in eq.(14) and 

eq.(15) we obtained the following equation[23]. 

r d(f) h 1 2 d(pi 1 - 2 d 2 (p h d 2 (pi T 

Ph{ — ^- + ] - ^[~^ + 2^') J + ~ p ^ = ^IP" ~ ^'a^ -ft + fh-fi (16) 

at the interface t) = r)(x,t). 

Substituting for 0^,0; from eq.(4),eq.(5) and value of using eq.(10)-eq.(12) and equating coef- 
ficient of x 2 , we obtain after some straightforward but lengthy algebraic manipulation, the following 
time development equation for 6 : 

rffe/v^) Nfa,r) (6/v^) 2 + 2(r _ 1} 6(66-1) (1 _ ^2*!) 



£>(&,r) (66-1) v 7 £>(6,r) v 

rc h [(s + l)(l-12e 2 2 )+46(s-l)] (17) 



0(6, r) 
where 

Pfc PMO k 2 v h ,2 (ph~pi)g nQ , 

r = — ;v h (i) = — = — ;c h = - 71 =;k c = 18 

Pi Ph(i) Ph vkg T 

Dfa, r) = 12(1 - r)il + 4(1 - r)£ 2 + (r + 1); 

7V(6, r) = 36(1 - r)il + 12(4 + r)& + (7 - r) (19) 

Eq.(8) and eq.(9) together with eq.(17) governs the temporal development of the Rayleigh- Taylor 
instability. For Richtmyer-Meshkov instability, the gravity dependent term in eq.(17) vanishes i.e. 
g — and the equation for ^ becomes 
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d(&/VZg) N(&, r) (6/v%) 2 24(r - 1)£ 2 3 (66 - 1) k 



dr D(£ 2 ,r) (66-1) Dfa,r) k 2 

] -rc h [(s + l)(l-12e 2 2 )+46( S -l)] (20) 



0(6, r) 



Ill SUPPRESSION AND OSCILLATION OF ASYMPTOTIC 



BUBBLE VELOCITY 

The nondimensionalized time development plots of £1,62 and £3 for RTI bubble are shown in Figure 

I. As r — > 00, the asymptotic values of £ 2 and £ 3 for bubble are obtained by setting ^ = giving 

d(4|=) 

[(&) aS ymp\bubbie = -| and — Jp 2 - = yielding 

" vV^ + fu^Mi-li) + i**> (21) 

It is interesting to note that if k 2 < 3/cj?, the asymptotic velocity of the bubble caused by the 
rising of the lighter (lower) fluid also by pushing through the denser (upper) fluid is affected only by 
the viscous drag exerted by the later. This is clearly seen from eq.(21) as [(£ 3 ) a symp]bubbie depends only 
on the kinematic viscosity of the denser (upper) fluid (vh,). If k 2 = 3k 2 equilibrium is attained; but 
if k 2 > 3k 2 , this reverse the sign of the second term in eq.(17) leads to the emergence of oscillatory 
state (Figure 2). 

Similar results for temporal development of Richtmyer-Meshkov instability are shown in Figure 
3. The asymptotic velocity is 



un i - 2k2 " h + 

[\<,3)asymp\bubble o 



4 9 A h 2 

^-kg——-)tm 



\9 h "9(1 + A)k 2 3 + A \9 h "9(1 + A)k 2 

For bubble as well as for spike, [{^) aS ymp\bubbie -» + \j\v\k^ - kg g( 2 ^ A) ^ exponentially 

with time. In absence of surface tension the time dependence is qualitatively similar to that for 
linear theoretical result of Mikaelian[17] but with different Atwoood number and kinematic viscosity 
coefficient dependence. However, the weakly nonlinear theoretical results of Carles and Popinet[20] 
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2 c h(i)kc^~A^ the asymptotic velocity of the bubble as well as the perturbed surface elevation oscillates. 



This is represented in Figure 4. 

In absence of surface tension both RTI and RMI are characterized by linear inviscid growth rate 
which increase with increasing wave number k. The dissipative effect due to viscosity also increases 
with increasing k and suppress the growth rate. A graphical representation of the wave number 
dependence of the nonlinear growth rate of the hight of the tip of the bubble is shown in Figure 
5. The nature of k— dependence is qualitatively similar to that in linear case [21] except that for 
t — > oo where the growth rate tends to a saturation value for all k in the nonlinear case. For bubble 



Thus for RTI the growth rate and perturbed interface are oscillating due to surface tension while 
for RMI oscillation depends on the relative strength of surface tension and viscous drag. 



To study the time evolution of spikes we adopt a procedure different from the usual Goncharov 
transformation[13]. We cast time evolution eq.(17) in the form given bellow. 



the asymptotic growth rate for RTI is maximum at [k max }bubbie 




IV TIME DEVELOPMENT OF SPIKE 



d(6/v^) 




dr 



\(t2 ~ l)rc h [(s + 1)(1 - 12^1) + 46(s - iMs/^fkg)} 



(23) 
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by using the following expression for N(£ 2 ,r) and D(£ 2 ,r): 
N(£ 2 ,r) = -36(r - 1)(6 - - «_) 

D(6,r) = -12(r - 1)(6 - - £-) (24) 

(r + 4) ± Vl6r + 4 
where ck± = — 



6(r- 1) 



-1 ± 



4r+2 



± = (25) 

Clearly a + > a_ and /3 + > /?_ and also a + > j3 + , a + > | 

First we consider the temporal development of spike as a result of RTI in an inviscid fluid with 
absence of surface tension, i.e., we put Ch = and = in eq.(17). We start with an initial value 
of 6 = 6o >«+>/?+ and |, and £ 30 < 0. 

Eq.(9) shows that ^ > while eq.(23) shows that d ^ 3 ^V / ^) < q ( as ^he square bracketed term 
on the RHS of the latter equation is positive) for all r when one starts from such initial values. Thus 
£2(7") is a monotonically increasing and £3/ \fkg is a monotonically decreasing function of time r. Now 
the curvature at the tip of the spike, i.e., x = is ^ = §^/[l + (§§ ) 2 ]^ = ^(t) = ^(t)- Thus the 
curvature of the spike increases with time while the acceleration of the tip of spike '^M^^l (directed 
downward) tends to a constant value as may be interred from eq.(23) when viscosity and surface 
tension are neglected. Thus the spike appears to fall continuously and simultaneously gets sharpened. 
This result is quite different from the earlier results obtained by Goncharov's transformation which 
concludes that the spike velocity tends asymptotically to a constant value. Our result is in conformity 
with expected spike behavior and is in agreement with some earlier theoretical results obtained by a 
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different approximate method [24] [25]. The same qualitative spike behavior is exhibited in presence 
of viscosity but with much reduced speed of fall of the tip of the spike. This is demonstrated in 
Figure 6 and Figure 7 obtained from the results derived from numerical integration of eq.(8), eq.(9) 
and eq.(23). Figure 6 shows that the spike speed decreases as the coefficient of viscosity increases. 
Moreover, for Ch = (inviscid fluid) the spike velocity is seen to vary linearly with time (i.e., close to 
free fall velocity) so that the displacement of the tip of the spike ~ \igt 2 where 7 is a dimensionless 
constant and value of 7 close to unity as A — > 1. 

Rayleigh- Taylor instability is driven by gravity g while Richtmyer-Meshkov instability is switched 
on by the impingement of a shock which impulsively changes the normal velocity by the amount 

= v a fter—Vbefore- Thus Richtmyer-Meshkov instability is driven by the instantaneous acceleration 
Av5(t). This has the consequence that the dynamical variables are to be nondimensionalized using 
normalization in terms of (kAv) for RMI instead of \fkg~ for RTI. Hence, in RMI equations -^=, Ch 
and r are replaced by 

W 9 = (k^r Ch = Jk&ry T = ( } (26) 

With replacements as gravity by eq.(26), the temporal development of RMI spike growth are 
obtained from numerical integration of eq.(8), eq.(9) and eq.(23) when the gravity g induced second 
term in the square bracket on the RHS of the last mentioned equation is to be deleted. The absence of 
gravity induced acceleration keeps close to its initial value when viscosity is neglected (c^ = 0) 
but tends to vanish when Ch 7^ 0. 

The lowering of the value of ^f^y in magnitude will according to eq.(9) reduce the rate of growth 
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of the curvature £ 2 with respect to r (nondimensionalized time). Thus, in contrast to the free fall 
of RTI spike, the RMI spike tip descends at an almost constant rate. The curvature also increases 
slowly, i.e, the spike sharpens at a much reduced rate as it falls (as compared to the RTI spike). All 
these features are shown in figure 8 and figure 9. The growth rate are seen to be further reduced due 
to the viscous drag. 

Finally, we note the difference between the RTI and RMI parabolic spike structures (figure 7 and 
figure 9) given by 

ky = krj (t) + {^)(kxf 

= Ut)+Ut){kxf (27) 

Let us first consider the inviscid case. The discussions in the foregoing paragraphs indicate that 

£i(t) ~ -\t 2 for RTI while ^(r) (constant) r for RMI (28) 

Further from eq. (8) and eq. (9) which gives (£ 2 + |) = (£20 + \)exp[— 3(£i — £ 10 )] (£10 and £ 20 are initial 
values) we obtain 

3 

€2(7) ~ (constant) exp[-r 2 ] for RTI while £ 2 (r) ~ (constant) exp[ (constant) r] for RMI (29) 

Eq.(28) and eq.(29) when plugged in eq.(27), now clearly indicates that the RTI spike fall much 
faster and gets sharpened much more rapidly. In both cases, the effect of viscosity is to dampen the 
growth rates. This leads to difference in the RTI and RMI spike structures as seen in figure 7 and 
figure 9. 
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V SUMMARY 



Finally we briefly summarize the results: 

(i) For RTI the asymptotic velocity of the tip of the bubble is given by 

K6)c 



I asympl bubble r. To — „ 

V /|^^ 4 +§TTTA)^(l-^)+§-^ 2 



and for RMI, the asymptotic velocity of the bubble is given by 

[{£z)asymp\bubble = — + ' \v\k^ — ^fi , 9 (i+ yl )"|j COth [( ^ ^/§^& 4 — % 9(1+^4) fj)^]- 

(ii) In case of RTI, if k 2 < 3k 2 the asymptotic velocity of the bubble is affected only by the viscous 
drag of the upper fluid(Figure 1) and similar effect for RMI with k 2 < 2c 2 h k 2 c ^^- (Figure 3). 

(iii) For RTI, the growth rate and perturbed surface are oscillating if k 2 > 3k 2 , i.e the oscillation 
depends only on the surface tension of the perturbed interface (Figure 2) while for RMI the growth 
rate and perturbed surface are oscillating if k 2 > 2c\k 2 ^-^-, i.e the oscillation depends on the surface 
tension of the perturbed interface as well as the coefficient of viscosity (Figure 4). 

(iv) In inviscid fluid the RTI spike has no asymptotically terminal velocity. Rather the spike has 
a nearly free fall so that the depth of the tip of the spike below unperturbed surface of separation 
~ \igt 2 where 7 is a dimensionless constant and value of 7 close to unity as A — > 1; also the spike 
sharpens as it falls. For viscous fluid the velocity of fall gets reduced as the coefficient of viscosity 
increases and tends to a nearly terminal velocity for sufficiently large viscosity. Similar result holds 
for RMI. The results are demonstrated in Figure 8 and Figure 9. 
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Figure 1: Variation of £i,£2 and fiz/sfkg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(17) for RTI bubble with initial value £1 = 0.1, £1 = —0.1, 
^j^m = 0.1, r = 1.5, jg = 0.5 and c h = (Red), 0.1 (Black), 0.2(Green), 0.3(Blue). 
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Figure 2: Variation of £i,£2 and fiz/sfkg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(17) for RTI bubble with initial value £1 = 0.1, £1 = —0.1, 
^j^m = 0.1, r = 1.5, c h = 0.001 and jg= 0.5 (Red), 5(Black), lO(Green), 15(Blue). 
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Figure 3: Variation of £1,^2 and fiz/sfkg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(20) for RMI bubble with initial value £1 = 0.1, £1 = -0.1, 
£3/^ = o.l, r = 1.5, jg = 0.5 and c h = (Red), 0.1 (Black), 0.2(Green), 0.3(Blue). 
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Figure 4: Variation of £1,^2 and fiz/sfkg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(20) for RMI bubble with initial value £1 = 0.1, £1 = -0.1, 
£3/^ = o.l, r = 1.5, c h = 0.001 and jg= 0.5 (Red), 5(Black), lO(Green), 15(Blue). 
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Figure 5: Variation of normalize £3 with normalize 'k' for RTI bubble with Ch 
r = 1.5, j£= and r = 3(Black), 8(Red), 00 (Green). 
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Figure 6: Variation of £i,£2 and ^,3/y/kg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(17) for RTI spike with initial value £1 = —0.1, £1 = 1, ^/\fkg = 
-0.1, r = 5,s = 1/5, jg = and c h = 0.0 (Black), 0.1 (Blue), 0.2(Red). 
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Figure 7: Formation of RTI spike where r = 5,s = 1/5, p- = and Ch= 0.0, r = 
2(Black), c h = 0.0, r = 2.5(Blue), c h = 0.1, r = 2(Red), c h = 0.1, r = 2.5(Green). 
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Figure 8: Variation of £i,£2 and fiz/sfkg with r as obtained by the solution of the 
eq.(8),eq.(9) and eq.(20) for RMI spike with initial value £1 = —0.1, £1 = 1, ^/y/kg = 
-0.1, r = 5, s = 1/5, jg = and c fe = 0.0 (Black), O.l(Blue), 0.2(Red). 
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Figure 9: Formation of RMI spike where r = 5,s = 1/5, §2=0 and c/,= 0.0, r = 
2(Black), c h = 0.0, r = 2.5(Blue), c h = 0.1,r = 2(Red), c fc = 0.1, r = 2.5(Green). RTI 
and RMI spike structures are different because the RTI is given by the acceleration g 
while RMI spike induced by shock velocity discuss below eq.(29) 



